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Abstract 



In this paper we study the competition between the Kondo effect and RKKY 
interactions near the zero-temperature quantum critical point of an Ising-like 
metallic spin-glass. We consider the mean-field behaviour of various physical 
quantities. In the 'quantum- critical regime' non-analytic corrections to the 
Fermi liquid behaviour are found for the specific heat and uniform static 
susceptibility, while the resistivity and NMR relaxation rate have a non- Fermi 
liquid dependence on temperature. 
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I. INTRODUCTION 



The interplay between Kondo screening of localised spins by conduction electrons, and 
ordering of these spins due to the RKKY interaction is a central issue in heavy fermion 
physics. Recently, a class of systems has been studied |l],FJ] in which the ordering temperature 
is driven to zero as a function of concentration, and the paramagnetic metal displays non 
Fermi liquid (NFL) behaviour near this T = quantum critical point. The compound 
Yx-xU x Pdz is one of the best documented among these systems 0,0] • In this case, the low- 
temperature ordered state is reported to be a spin-glass for x > x c ~ 0.2 , while the system 
remains a paramagnet down to the lowest temperature studied for x < x c . It is still debated 
HH whether the NFL behaviour of the lo.8^o.2-P^3 system is a single-ion effect or results 
from the above competition and the proximity of the T = critical point. The aim of 
this paper is not to resolve this debate for this particular system, but to demonstrate that 
NFL behaviour is indeed a generic feature of the vicinity of a T = paramagnetic metal to 
metallic spin-glass transition. This will be shown by solving specific models at mean-field 
level. 



II. MODELS 

The models that we shall study are mean-field versions of the Kondo lattice, with an 
additional quenched randomness on the exchange interactions between localised spins. We 
consider localised spins Si on a fully connected lattice of N sites % = 1, • • • , N. These spins 
interact with a bath of conduction electrons. In the model that we shall consider first, a 
major simplification will be made: the conduction electron bath will be assumed to consist 
of independent "reservoirs" of electrons, with no communication between the reservoirs at 
different sites. The effect of releasing this simplifying assumption will be discussed in Sec.[TV| 
at the end of this paper. The conduction electrons will be denoted by cj^,, where a =f, J, is a 
spin index, k labels the conduction band orbitals, and the site index (i) labels the reservoir 
associated with site i. The hamiltonian of the model reads: 

Hi = E e * c M + J '< E & ■ - E W ■ s- (i) 

k,a,i i i<j 

In this expression, s(i) = J2 a /3 J2k,k' \ c ka &a/3 c k>B * s the conduction electron spin-density 
at site i, and Jk the strength of the Kondo coupling between the localised spins and the 
conduction electrons (taken to be antiferromagnetic). Besides this coupling, the spins have a 
direct interaction between one another: the J^'s are quenched independent random variables 
with a distribution P{Jij) ~ exp(— JfjANJ 2 ). A further simplification of our model is that 
only the Ising part of the exchange interaction has been included. For Jk = 0, the model 
reduces to the Sherrington-Kirkpatrick model of a classical Ising spin-glass || with a freezing 
transition at T c = J. In contrast, for J = 0, we have a system of independent localised spins, 
each one being quenched by the Kondo effect with its conduction electron reservoir, and the 
system has no long-range order down to T = 0. We are interested in the intermediate 
behaviour where the spin-glass freezing due to the random exchange competes with the 
local Kondo effect. 
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Because the model is fully connected, it can be reduced to a single-site problem after 
taking the (quenched) average over the realizations of the random couplings. In order to 
describe also the spin-glass phase, we may introduce replicated variables labeled by indices 
a,b = 1, • • • , n. Using standard techniques [|7|-§|], the single-site effective action is found to 
be (with obvious notations): 

Seff = ~ SS dr SS dr> L, a c+ a (r)g o l (r - r')c Q , a (r') + J K dr £ a S a (r) ■ ^(r) 

- J 2 Ji dr Ji dr' Ea, b S z a (r)D ab (r - r')^(r') (2) 

In this equation, Qq is simply the bare on-site propagator for the conduction electrons: 
Qq = J2k l/(^n — ejfc)j which can be taken to have the characteristic form associated with a flat 
band without loss of generality: £ (?c<j n ) _1 = iTsgn(u) n ). The non-local spin interaction D ab 
must satisfy the following self-consistency condition in the thermodynamic limit iV — > oo: 

D ab (r - t') =< TSI{t)SI{t') > Seff (3) 

In the following, we shall mainly be concerned with the paramagnetic phase in which replica 
symmetry holds, so that D ab (r) is non-zero only for a = b. In the spin-glass phase, D ab 
is non-zero (but r-independent) for a ^ b, and the Edwards-Anderson order parameter is 
given by q EA = D m (r -> oo). 

The problem defined by Eqs.(fJfJ) is still too complicated to be solved exactly, even in 
the paramagnetic phase. However, it can be related to a solvable problem, which can be 
shown to have the same qualitative phase diagram and the same low-frequency and low- 
temperature universal properties. To arrive at this solvable model, we have to go through 
two additional steps. The first one is to "integrate out" conduction electrons in S e ff, so 
that an action involving only spin degrees of freedom is obtained. This cannot be done 
exactly because of the Kondo interaction, but can be done asymptotically at low-energies 



by following the classic Anderson- Yuval-Hamann approach to the Kondo problem [10]. This 
approach consists in separating the Kondo term into an Ising part S z s z and a spin-flip part, 
and performing an expansion in the spin-flip term to all orders. The result of this procedure 
is a mapping of the Kondo part of S e ff onto an action involving an effective interaction for 
the Ising components of the spins. This effective interaction is retarded and decays as 1/t 2 
for large (imaginary time) separations r. Thus, low-energy properties of the model can be 
studied by replacing S e ff with: 



S' eff = fdr f dr'Y.S z a {r)[K{r - r')6 ab - J 2 D ab {r - r')\S z h {r') (4) 



in which K(t) ~ 1/t 2 at large r, i.e K(iuj n ) ~ K,u n sgn(u) n ) (« is a dimensionless parameter 
depending on the anisotropy of the Kondo coupling). This behaviour must be cutoff at 
high frequency uj > A (i.e at short time separation). Because the decay K(r) ~ 1/t 2 only 
holds for times up to 1/Tk, where Tk is the 'bare' Kondo temperature in the absence of the 
exchange J, the cutoff must be chosen as: A ~ Tk- This slow decay of K(t) stems from 
the metallic nature of the spin-glass problem under consideration. In the insulating case 



considered in [11 , 12| , the decay is slower with K(iuj n ) ~ uj. 
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III. SOLUTION AND PHYSICAL PROPERTIES 



A. Solvable M = oo limit 



In order to analyze the effective action (£|), we can follow Ye, Sachdev and Read []T2| and 
generalize the model from Ising spins to M-components quantum rotors n ifl , \x = 1, • • • , M, 
with nf = 1 at each site. The Ising case formally corresponds to M — 1, but the M — >• oo 
limit will be first solved exactly Enforcing the local constraint by Lagrange multipliers Xj, 



the model is solved in this limit by a saddle-point method ||12|| , with iXj = A uniform at the 



saddle point. The free energy of the model reads, in the paramagnetic (replica symmetric) 
phase: 

im = x + h^ lYl[x + K{lUn) ~ j2D{lUn)] (5) 

where the spin-correlator D(iuj n ) is given by: 

1 



Dfan) = ^pi* + %n) - V( A + K (^n)) 2 ~ 4J 2 } (6) 

The Lagrange multiplier A is determined by the constraint equation (equivalent to setting 
dF/d\ = 0): 

-^XX^n) = l (7) 



B. Connection with z = 2, d = 3 quantum critical phenomena 

In order to analyze the phase diagram and critical behaviour resulting from these equa- 
tions, it is very useful to put them in a different form, which will reveal a connection with 
a different problem already analysed in the literature. Indeed, instead of using a mapping 
onto the single-site action (§), we could have solved the M = oo rotor model directly for 
the fully-connected lattice. In this approach, the spin correlation function D(iu n ) is given 
by the on-site component of the inverse of the random matrix: (A + K(iu n ))Sij — Jij. In 
the thermodynamic limit N = oo, the eigenvalues of the matrix have a semi-circular 
distribution given by: 

1 



Pj(x) = v / 4J 2 -x 2 ^(2J - \x\) (8) 

Hence, the free-energy and spin propagator can be written in the equivalent form: 
F r+ 2J 



MN 



/ + ZJ 
dxpj(x) ln[A + K(iu n ) - x] (9) 

2J A + Kiiujn) — X 



D{lU ri ) = l +2J dx- (10) 
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Independently of the low-frequency behaviour of K(iu n ) (metallic or insulating), it is clear 
that only the solutions with A > 2J are admissible []12| (so that no pole is encountered 
in the integration over x). The critical boundary with the spin-glass phase is signalled by 
A(T, J) reaching the value 2 J (the condition A(T, J) = 2 J holds throughout the spin-glass 
phase). Thus, the combination A = A — 2 J plays the role of the important low-energy scale 
(which vanishes in the spin-glass phase). Universal low- frequency properties are found when 
this energy scale is small. These universal properties only depend on the low-frequency 
behaviour of K{iu) n ) and on the fact that the spectral density pj(x) has a square-root 
behaviour near its upper band edge +2 J. Formally, this behaviour is identical to that of 
fermions with a kinetic energy proportional to k 2 in d = 3 dimensions: we could as well set 
K(iu n ) + A — x = K{iuj n ) + A + k 2 and replace the integration over x in the above equations 
by an integration over d 3 k (with some upper cutoff). In this analogy, the scaling of frequency 
with respect to k is obtained from K(iu n ) ~ u 2 ~ k 2 in the insulating case of Ref. |12| 



while K(iu; n ) ~ \u n \ ~ k 2 in the metallic case. Therefore, we conclude that there is a formal 
equivalence with quantum critical phenomena [|n|,|14[] in the d = 3, z = 1 universality class 
for the insulating case, and in the d — 3, z — 2 universality class for the metallic case of 
interest here. The insulating case considered in |L2| corresponds to a (quantum) Landau- 



Ginzburg model at its upper critical dimension d + z = 4 (with 4 a marginal perturbation), 
while the metallic case amounts to look at this model above its upper critical dimension 
(with <f) A a dangerously irrelevant perturbation). For these reasons, many results that will 



be derived below for the metallic case are formally identical to those of Millis |14| for the 
d = 3, z = 2 case. 

In order to derive these results, it is most useful to convert the Matsubara sums in the 
above equations into real-frequency integrations. Using K(lo + i0 + ) ~ iu sgn{uj) (with an 
upper cutoff A), we can write the free energy as: 



NM 

and the constraint equation reads: 



A + / — coth — / dxpj(x) tan (11) 

Jo n 2 J-2J A — x 



A duj X "(oj)coth^ = 1 (12) 
o 2 



with: 

7T J-2J 7T[UJ Z + (A — xy\ 



1 f+ 2J 10 
X'{uj) = lmD(u + iO + ) = sgn(u) / dxpj{x) — (13) 

7T J-2J 7T \UJ Z + A — X r 



C. Phase diagram 

The phase diagram resulting from Eqs. ( [I2]JT3] ) is depicted schematically in Fig.|I| It is 
qualitatively similar to the insulating case |12|, but the equations for the critical boundary 
and the various crossover lines are affected by the different low-frequency behaviour of 
K(iuj n ). Some calculations are detailed in the Appendix. For very large J, the spin-glass 
transition temperature is at T c = J, the classical value ||. T c decreases upon increasing 
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quantum fluctuations (i.e increasing Jk) and eventually vanishes for J smaller than a critical 
value: J c ~ An 1 / 3 ~ T^k 1 / 3 . Near zero temperature, the phase boundary is such that: 
1 - J C {T)/J C ~ (T/A) 3 / 2 ~ {T/T K f/ 2 (to be contrasted with ~ T 2 in the insulating case 

m 



Next, we discuss the various crossover regimes |TJ]] near the T = quantum critical 
point at J = J c (cf. Fig.p]). Close to this point, and for low frequency and small A 
(but uj/A arbitrary), the imaginary part of the local dynamical susceptibility x"^) = 
— l/irImD(u + i0 + ) takes the scaling form: 



° sgn{u)^\f(^ (14) 



Jl' 2 

The scaling function f(x) is easily obtained from (|(J) as: 

f(x) = (x/{i + VT+^)) 1/2 (15) 

and behaves as f(x) ~ 1 for x —>■ oo and f(x) ~ y/x for x — > 0. Hence, spin excitations are 
gapless in all regimes, with x"^) ~ ^/(^c^v^A) f° r w << A an d ~ sgn{oj)\/uj '/ ' J^ 2 

for u; >> A. The former, linear behaviour is characteristic of local spin correlations in a 
Fermi liquid with a low-energy scale A, while the latter (which holds down to uj = at 
the T = critical point) deviates strongly from Fermi-liquid theory. C is a non-universal 
constant, and the low-energy scale A(T, J) = X — 2J has different behaviour in different 
regions of the phase diagram (Fig.|l]). 

A = inside the spin-glass phase, in which x"^) ~ sgn{uj)\/uj with no characteristic 
scale. 

Raising temperature at J = J c , one enters the quantum critical (QC) regime, in which 
the physics is dominated by the T = quantum critical point. The energy scale A = A — 2 J 
is set entirely by temperature in this regime and is found to be A ~ T 3 / 2 /^T c ~ T 3 / 2 /^T^. 
One would have naively expected A ~ T but we find a violation of this scaling in the QC 
regime of this model. Following the connection explained above, this is an effect of the 4 
term in the quantum Landau- Ginzburg model being a dangerously irrelevant perturbation 
for d + z = 5 > 4. 

Decreasing J (or increasing T%) from the T = critical point, at low enough temperature, 
one enters the quantum disordered (QD) regime in which A ~ J c — J ~ T# — J (to dominant 
order). The crossover between the QC and QD regimes occurs at J = J*(T) obtained by 
matching the two behaviour of A given above, with the result: 1 — J*(T)/J C ~ (T/J c ) 3 ^ 2 ~ 
(T /Tk) 3 ^ 2 ■ As shown below, two distinct regions corresponding to T ^> A and T < A 
must actually be distinguished within the QD regime, in which the physical quantities have 
quite different low-temperature behaviour. In the low-temperature region of the QD regime 
(denoted QD 2 on FigfJ), the physics is that of a metal showing Fermi-liquid behaviour below 
the coherence scale (or effective Kondo scale) A, which can be very small because of the 
competition between the Kondo effect and the freezing of the local moments. However, this 
interpretation has to be handled with some care, since this scale enters the various physical 
quantities in quite different manners, as detailed below. 

Finally, near the phase boundary, there is a classical regime jl4j for \J — J C (T)\ << T 2 /J c 
dominated by purely classical fluctuations and in which A ~ J c ( J — J C (T)) 2 /T 2 . 
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D. Specific heat, Susceptibility and NMR relaxation rate 

We now investigate the low-temperature behaviour of the specific heat in the various 
regimes described above. Some details of the calculation are provided in the Appendix, but 
the results could also be directly read off from those in [fL4}| , given the equivalence with the 
d = 3,z = 2 problem. In the QC regime, the specific-heat coefficient 7 = C/T is found to 
behave as: 

7 = 7o--^v / T + --- (16) 

Jc 

Interestingly, the non-Fermi liquid nature of the quantum critical point results in a non- 
analytic correction to the low-temperature specific heat in this regime. Very close to the 
critical boundary with the spin-glass phase, there is an additional contribution ||14j| : BT/[J— 
J C (T) + T 3 / 2 ] 1 / 2 (~ T 1 / 4 for J ~ J C (T)), which becomes rapidly negligible however as one 
moves away from the critical boundary. 

The non-analytic \/T behaviours continues to hold within the QD regime as long as T ^> 
A, i.e T ^> J c — J. This defines an additional subdivision of the QD regime, corresponding 
to region QD : in Fig.|l|. In the low-temperature (QD 2 ) region of the QD regime, (|16"D is 
replaced by: 

C D T 2 

7 = 7o - -3/2 VA - ^ + • • • (17) 

with A ~ J c — J and 7 oc 1/J C remains finite as J — ► J c . 

Next, we discuss the low-temperature behaviour of the uniform spin susceptibility \- 
This quantity is assumed to be measured by applying a uniform field such that the Zeeman 
energy \xH is smaller than the bare Kondo scale Tk- Then, the Kondo effect still takes place 
and the low-frequency behaviour of the resulting propagator K{iui n ) is essentially unaffected. 
We can thus simply introduce a uniform magnetic field HJ2i Sf in the effective Ising model. 
The susceptibility is given by the sum over all pairs of sites of the spin-spin correlation 
function taken at u n = + , namely: namely: 

X = ^[A + AT(iO+)- J w ]^ (18) 

u 

In this expression, the overbar denotes an averaging over disorder. The inverse of the random 
matrix present in this expression can be evaluated by expanding in powers of Jm- Because 
these couplings have random signs and zero- mean, only closed paths (with i = j) give a non- 
zero contribution on the fully-connected (or Bethe) lattice. Hence, the uniform susceptibility 
behaves in an identical manner to the local spin susceptibility: xioc — Jc \ drD(r). Note that 
this relies crucially on the fact that we are dealing with a random system, and would not 
apply to a a uniform antiferromagnet for example. Hence the formal analogy between the 
present problem and the z = 2, d = 3 quantum antiferromagnet does not apply to the 
calculation of the susceptibility, xioc (and thus x) is easily obtained by setting uj n + in 
Eq.(|) and taking the real part. This yields 2J 2 ^ oc = const. - 2(JA(T, J)) 1/2 + A(T, J). 
In the QC regime, we must set A ~ T 3//2 , so that: 



7 



x(T) = X o(l-a(?-f '* + ■■■) (19) 

" c 

with Xo ~ l/^o an d c a numerical constant. A non-analytic correction departing from 
standard Fermi-liquid theory is again found. In the high-temperature part of the QD regime 
(QD 1 region), we have (cf. Appendix): A(T) = A + T 3 / 2 , so that: 

X (T) = Xo - ai -^-a 2 j^= + --- (20) 

In the low-temperature (Kondo) region of the QD regime, we have ( cf Appendix) : A(T) = 
A + T 2 /v / A, so that: 



V A T 2 

x{ T )=Xo-ai— a2 j2^ + '" ( 21 ) 

In both ( PU|) and ([21]), A ~ J c — J and we emphasize that %o ~ l/Jc is non-singular as the 
critical point is reached. 

Finally, we consider the NMR relaxation rate, which is directly related to the u = 
behaviour of the local dynamical susceptibility through: 

1 -*" M U„ (22) 



Hence, this quantity always feels the linear regime x"^) ~ cj/v^A in the scaling form (|Tl|). 
In the QC regime, this yields a temperature dependence of the relaxation rate which differs 
from the usual Korringa behaviour (1/T\T ~ const.) found in a Fermi liquid: 

1 1 (23) 



In contrast, in both the QD X and the QD 2 ('Kondo') regime, the Korringa law is obeyed, 
but with an enhanced rate: 

1 1 1 (24) 



We emphasize that in the same regime, the T = uniform spin susceptibility is not critically 
enhanced (the low-energy scale A only enters subleading corrections). 



E. Ising (M = 1) case 

Before leaving model (|l|), we would like to show that the Ising case M — 1 actually has 
the same universal low-frequency behaviour than the rotor model in the M — > oo limit that 
we have analysed in detail. In order to show this, one possibility is to adapt the reasoning of 
Refs. [0,0 to the present case. Specifically, we can define the spin irreducible self-energy 
associated with (j|) by: Diiujn) -1 = K{iuj n ) — J 2 D — n(iu; n ). For M = oo, H(iu> n ) reduces 
to a constant II = A. To order 1/M, the first contribution to ImU is from decays into 



S 



three spin waves. With D{t) ~ 1/t 3 / 2 at the critical point, this leads to ImH ~ ^J/ 2 - 
No non-analyticity with a weaker power of frequency is induced to any order in 1/M, and 
hence the low-frequency behaviour found above at the critical point is unchanged. Another 
line of reasoning leading to the same conclusion is to use the formal equivalence with the 
d = 3, z = 2 universality class. Going from the M = oo to the M = 1 case just changes the 
specific coefficients of the various terms of the (quantum) Landau-Ginzburg model, but the 
equivalence holds for all M. 



IV. MODEL WITH A SINGLE CONDUCTION ELECTRON FLUID AND 

TRANSPORT PROPERTIES 

Finally, we shall release the assumption of independent Kondo baths at each site, and 
consider a model with a single species of conduction electrons. We shall consider a lattice 
of connectivity z, and study the hamiltonian: 

H 2 =-]T UjC+Cjv + J K Y, Si ■ s(i) - J2 J nSi ■ SJ (25) 

ija i <ij> 

The random exchange couplings are distributed according to a gaussian distribution as 
above, which is now normalized such that Jfj = J 2 /z. The limit of large lattice connectivity 
z — > oo will be considered below. In this limit, the various Green's functions satisfy self- 
consistent dynamical mean- field equations H|J. These equations reduce the model to the 
solution of a single-site effective action, which has again precisely the form in Eq.Q. However, 
there is now a self-consistency condition on both D{iuj n ) (given by Eq.|| above) and Qo(iu n ), 
which is not known explicitly in contrast to the case above. This self-consistency condition 
depends on the specific lattice and hopping term Uj. More precisely, if D(e) stands for the 
non-interacting density of states of the lattice under consideration, the effective propagator 
Go must be such that the following self-consistency equation holds |||9|: 

/"+ 00 Die) 
G c ^ n ) = / ale—— { -L— (26) 

In this equation, G c stands for the local conduction electron Green's function G c = — < 
Tcc + >seff, and S c for the self-energy S c = Q ~ 1 —G~ 1 . Both should be viewed as functionals 
of Qo{iuj n ) (and D(iu) n )). 

A first possibility is to consider a model with a specific form of long-range hopping 
(described in 0), such that the non- interacting conduction electron density of states is a 
Lorentzian D(e) = l/vr(e 2 + T 2 ). In this case, the Hilbert transform in (|26| ) yields: G c = 
[iu n + n — S c + iTsgn{oj n )}~ 1 , so that S c disappears altogether from the self-consistency 
equation and Go is actually known explicitly as before: Qq 1 = ito n + iTsgn(u n ). Hence, 
exactly the same equations as in the above model (P with independent Kondo baths are 
found for this model with long-range hopping, and no non-trivial feedback of the conduction 
electron dynamics into the spin dynamics is possible. 

The situation is different for a model with short-range hopping of the conduction elec- 
trons. For definiteness (but without loss of generality), we may consider the z = oo Bethe 
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lattice, with nearest-neighbour hopping normalised according to ty = t/yfz. (This corre- 
sponds to a semi- circular d.o.s with a half- width 2t). In this case, (|26f) takes the simpler 
form: 

G^iioJn) =iuj n + fi- t 2 G c (iuj n ) (27) 

G c has to be determined from the solution of S e ff itself, so that the problem involves a 
self-consistency condition on both D{iu n ) and Go(iu n ). The main question is whether this 
"feedback" of the non-trivial spin dynamics into Qq can change the low-frequency behaviour 
of x"(u;) at the T = critical point. We shall give an argument that this is not the case, 
and that (0) still holds. 

Let us imagine that the coupled problem (p|,|3|j27|) is solved iteratively, starting from 
a Qo which has the same long-time behaviour than in the model without feedback, namely 
Go( T ) ~ 1/r. Then, the arguments above yield D{t) ~ 1/r 3 / 2 at the critical point. Inserting 
this into the effective action (Q), we have to compute the conduction electron Green's function 
G c and feed it back into the self-consistency condition ( P7D to see how Q is affected. In 
order to find the behaviour of G c , it is convenient to use a representation of the localised 
spins by pseudo-fermions f a , such that: S = f^a aa / f a /2. This amounts to "undoing" the 
Schrieffer- Wolff transformation and representing the original Kondo lattice as a periodic 
Anderson model with a very large U in the local moment limit ej ~ —U/2. We shall 
not attempt here to find the actual behaviour of the /-electrons Green's functions Gf(r), 
but it is easily seen that it cannot decay faster than l/r 3//4 . Indeed, D(r) always contains 
a term of the form Gf{r)Gf{— r) (supplemented by vertex corrections), so that a decay 
slower than l/r 3//4 is inconsistent with D{r) ~ 1/r 3 / 2 . Furthermore, the slowest possible 
decay of G/(r) is the Fermi liquid form G/(r) ~ 1/r. Hence, we conclude that Gf(r) ~ 
l/r e , with the exponent 9 such that 3/4 < 9 < 1. Now, the t-matrix associated with (Q) 
is proportional to Gf(iu n ), so that the conduction electron Green's function is given by: 
G c {iu n ) = Q (iuj n ) + V 2 Q (iu n ) 2 G f(iuj n ) . Hence G c {uj) behaves as u 8 ^ 1 at low-frequency. 
Inserting this into the self-consistency equation fl27l) in order to see how Q is affected at 
the next step of the iteration, we see that Qo{r) behaves as l/r 2 ~ 8 for large r. Because of 
the constraint above, this exponent satisfies the bounds 1 < 2 — 9 < 1 + 1/4, so that Go 
cannot decay faster than 1/r at the next iteration. Hence, the associated spectral density 
ImGo{uo + i0 + ) cannot diverge at u = 0: it is either finite or vanishes. For large enough 



Kondo coupling, this yields a standard Kondo effect [13 , and hence an effective interaction 
between Ising spins such that K{r) ~ 1/r 2 as before, so that the behaviour of ~x!'{uj) at 
the critical point remains unchanged at any step of the iterative solution of the coupled 
equations. 

A perturbative argument can actually be given that the dominant behaviour G c (r) ~ 1/r 
of the conduction electron Green's function is not affected at the T = critical point, i.e 
that 9 = 1. Indeed, if we treat the residual coupling between the conduction electrons and 
localised spins in second- order perturbation theory, we obtain a contribution to the self- 
energy: E c (r) oc J|-D(r)^o(' r ) ~ 1/t 5 ^ 2 , so that S c (o;) behaves as cu 3 ^ 2 at T = 0. Hence, this 
scattering is unable to modify the dominant term in the long-time behaviour of G c (r) ~ 1/r. 

It has however important consequences for the transport properties near the T = 
critical point, as we shall now show. We first perform a more precise evaluation of the 
finite-temperature scattering rate to order J^, which takes the form: 
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ImS c (^ + ,0+) x 4 Po (0) du X »(u) + ^ + J (28) 

Inserting the scaling form fll4|) , this leads to the following low-frequency and low-temperature 
dependence, in the quantum critical regime: 

Im£ c (cu + *0+) oc ^ po (0)[o// 2 + T 3 / 2 ] (29) 

Jc 

while in both regions of the QD regime: 

Im£ c (^ + *0 + ) x j& po(0)[-^ + ^] (30) 

Since we are dealing with a model on a lattice with infinite connectivity, the vertex correc- 
tions to the conductivity vanish |IB |, and the (ic-conductivity is simply given by: 

a dc oc / deD(e) / du - — ■ — — ^—77 v , 2 31 

Hence, the above calculation of the scattering rate leads to the following non-Fermi liquid 
temperature dependence of the resistivity in the QC regime : 

5p(T) ~ T 3/2 (32) 

while in both regions of the QD regime: 

T 2 

Sp(T) ~ -= (33) 



with A ~ J c — J. Hence a Fermi liquid behaviour 5p = AT 2 is recovered in the QD regime, 
but with a critically enhanced rate A ~ 1/ y/A. We emphasize that in the same regime, the 
specific heat coefficient 7 is not critically enhanced. This distinguishes the QD Fermi-liquid 
regime from conventional heavy-fermion behaviour in which both A and 7 are large for small 
T K , with p| A oc 7 2 . 



V. CONCLUSION 

In this paper, we have studied Kondo lattice models with a quenched random exchange 
between localised spins. The mean-field phase diagram has been investigated (Fig.|l|) and 
found to display several different regimes near the quantum critical point associated with the 
T = spin-glass transition. In the 'quantum critical' regime, the specific heat coefficient and 



susceptibility display non-analytic corrections to Fermi liquid behaviour given by ( p^|T9"D , 
while the NMR relaxation rate (|23| ) and resistivity (|32]) have a non- Fermi liquid temperature 
dependence. In this regime, the important low-energy scale violates naive scaling and varies 
as a power of temperature (A ~ T 3 / 2 at mean-field level). In the low-temperature part of 
the quantum-disordered region ('Kondo regime'), Fermi-liquid behaviour is recovered, but 
the NMR and scattering rate are critically enhanced as the transition is reached, while 7 
and x are n °t- 
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These results may have qualitative relevance for Yi^ x U x Pd 3 and related systems 
since they indicate that non-Fermi liquid behaviour is a rather generic feature associated 
with a T = spin-glass transition in a metallic system. However, the reported experimental 
behaviour (7 ~ — InT, X = Xo ~ VT, 5p ~ T) is not in good agreement with our mean-field 
results. This raises theoretical questions associated with the fluctuations beyond mean field, 
and also experimental questions concerning the actual investigation of the critical scaling 
regime. 
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APPENDIX A: 

In this Appendix, we provide some details on the analysis of the various crossovers and 
on the calculation of the specific heat. The starting point are the Eqs.(fn,l2,I5D, in which 



we set: A = A — 2 J and change variables in the integrations over x by setting x = 2 J — e. 

We shall deal first with the constraint, Eq. QTi2"|) , which reads (replacing pj(x) by its 
square- root form near the upper band edge): 

f A dujf 4J de^Buj uj 
1= — — ^v^coth^- Al 
Jo 7r Jo 7iJ 6 ' z 2 uj z + (e + A) z 

The integral over uj needs to be cutoff at A ~ Tk, while the integration over e is ultra-violet 
convergent and its upper limit could as well be set to +oo. We shall focus on the regions 
where A < T (which includes the QC regime and the upper part of the QD regime). Under 
this assumption, an expansion in A can be performed to yield: 

A duj f 4J de ,- , 3u ( uj . 2uje ^ /A 9n\ , , „ N 

1 vWh^- - - A : ■ ^ + 0(A 2 ) (A2) 



o 7T jo Ti.Pl 2 v 2 \u 2 + e 2 (u 2 + e , 

In the integrations over e, we set e = uoj so that the uj dependance becomes apparent. Then, 
we use coth/?c<j/2 = 1 + 2/(e w / T — 1) and expand the integrals involving the last term at low 
temperature. The constraint equation then takes the form: 

i = M j) 3/2 + M j) 3/2 - j ( 5 i(j) 1/2 + ^(j) 172 ) + °( A2 ) ( A3 ) 

In this expression, the Aj's and B^s are purely numerical, positive constants (independant 
of the cutoff and of J). 

The location of the quantum critical point is readily obtained by setting A = T = 0, 

2/3 

yielding: J c = A{ A oc T%. When T is increased above this point, the behaviour of A 
in the QC regime is found by cancelling the dominant terms to next order, leading to: 
A ~ Ail B{F Z I 2 /a/A. The above expansion is valid as long as A < T, so that we can also 
use it for T 3 / 2 <C ACT, corresponding to the upper part of the QD region. In this case, we 
have to expand in J = J c — 5 J and the behaviour A ~ 5 J + 0(T 3//2 ) is obtained. Also, the 
shape of the critical boundary at low temperature is found by expanding in 5 J = J C (T) — J c , 
with A set to 0. This yields 5 J oc T 3 / 2 /v% 

Next, we give some indications on the low-temperature expansion of the free-energy in 
the QC regime. We rewrite Eq . (|TT|) for the free-energy per site / = F/NM under the form: 
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/ = const. + / — coth— / — - -Wetan — — (A4) 

In the QC regime, we should use the following scaling variables: uo = Too, e = Te, A = T 3 / 2 A. 
Thus it is clear that A is the smallest energy scale, and we can simply expand the above 
expression in powers of A. The important point here is that the linear term vanishes because 
of the constraint equation above. Hence, this expansion reads: 

f A doo , Boo f 4J de ^ ( ,u a9 ooe ^ /A ^\ 
/ = const. + / — coth^- / — — ^ tan' 1 - + A 2 + O A 3 A5 

Jo i 2 Jo 7rJ d / 2 y e (or + e 1 ) 1 J 

The coefficient of the A 2 term is the same as found above in the analysis of the constraint. 
It leads to a dependence of the form A 2 (const. + VT) ~ T 3 + T 5/2 . 

In order to find the temperature dependence of the first term, we make a low frequency 
expansion: 

L ^^tan-'-~ 7 + ( 7 P + 0(^) (A6) 

Hence, this term yields a contribution of order T 2 + T 5 / 2 to the free-energy. Overall, we find 
the low-temperature expansion: 

/ ~ const. + }\T 2 + f 2 T b ' 2 + ■■■ (A7) 

in which the coefficients /j's are non-singular and non-vanishing as the critical point J = J c 
is reached. The behaviour of C — Td 2 f/dT 2 ~ 70T — T 3 / 2 follows. 

Finally, we comment on the behaviour of A in the low-temperature QD 2 regime. In this 
regime, the above expansions are no longer valid since they assumed T ^> A. Concentrating 
e.g on the constraint equation, we see that the leading low-temperature correction is now 
controlled by the long-time behaviour D{t) ~ l/r 2 V / A, which holds for rA ^> 1. Using the 
Poisson summation formula, this yields: 

11 T 2 

- D{iuj n )\ J= j C! T=o - „ E D ( iuJ n) = a5J + 65 A + c^= + • • • (A8) 

P n P n VOJ 

with 5 J = J — J c , 5\ = A — 2J C . Cancelling the leading corrections, one obtains: A = 
j _ j j Tl 
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FIGURES 

FIG. 1. Schematic phase diagram at mean-field level, as a function of T/J and J c /J (with 
J c ~ A ~ Tk)- The plain line is the critical boundary with the spin-glass phase. All other lines 
are crossover lines, corresponding to the regimes described in the text. The hatched region is that 
of classical behaviour near the critical boundary. 
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